Suppression of core polarization in halo nuclei 
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. We present a microscopic study of halo nuclei, starting from the Paris and 

Bonn potentials and employing a two-frequency shell model approach. It is 
found that the core-polarization effect is dramatically suppressed in such nu- 

Ij . 

^ . clei. Consequently the effective interaction for halo nucleons is almost entirely 

given by the bare G-matrix alone, which presently can be evaluated with a 
high degree of accuracy. The experimental pairing energies between the two 
halo neutrons in ^He and ^^Li nuclei are satisfactorily reproduced by our cal- 
culation. It is suggested that the fundamental nucleon-nucleon interaction can 
be probed in a clearer and more direct way in halo nuclei than in ordinary 
nuclei. 

PACS numbers: 21.30. -hy,21.60.Cs, 21.90.-hf 
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Radioactive-beam nuclear physics has been progressing rapidly, and there is much current 
interest in studying halo nuclei There were four articles about halo nuclei in a recent 

issue of Physical Review C: Nazarewicz et al. dealt with the halo nuclei around the nucleus 
^^Ni, Hamamoto et al. 0] carried out a systematic investigation of the single particle and 
collective degrees of freedom in the drip-line nuclei, and an experimental study of heavy halo 
nuclei around = 82 was also presented It is remarkable that nuclei as exotic as ^^Ni, 
i.e., the mirror image of the double closed-shell nucleus ^t'a, are now being studied! 

The halo nuclei (or drip-line nuclei) may well play a central role in our understanding 
of the nuclear binding. Their typical structure is that of a tightly bound inner core with 
a few outer nucleons that are loosely attached to the core. Although these exotic nuclei 
are bound, their binary subsystems are not. For instance, the halo nucleus '^He {}^Li) is 
presumably made of a ^He i^Li) core surrounded by a two-neutron halo. As a whole ^He 
i^^Li) is bound, but its binary subsystems, i.e., ^He {^^Li) and the di-neutron are unbound. 
The pairing force between the valence nucleons is thus essential for the stability of the halo 
nuclei, and it is important to calculate it as accurately as we can. 

So far, the halo nuclei have been calculated using empirical effective interactions, tuned 
to stable nuclei. The inherent density dependence of Skyrme-type forces [^,^ provides a 
reasonable means of extrapolating to the lower density regimes characteristic of nuclei far 
from stability. Yet, quite recently Kuo et al. have suggested to study the drip-line nuclei 
from the first principles, i.e., from the elementary nucleon-nucleon (NN) force, such as the 
Paris 0] and Bonn p]] interactions. Halo nucleons are separated rather far from the other 
nucleons in the "core nucleus", and the interaction among them should be derivable from 
the free NN interaction with small medium corrections. 

Starting from a free NN interaction, a model-space effective interaction (K//) among the 
nucleons in the nuclear medium can be derived, using a G-matrix folded-diagram approach 



PJTO[]. The major difficulty in such a microscopic effective interaction theory has been 
the treatment of the core polarization effect (CPE), in particular the higher-order core 
polarization diagrams. In ordinary nuclei, the valence nucleons are close to the nuclear 



core, an example being the two sci-shell neutrons of residing adjacently to the 
core. Consequently, there is a strong valence-core coupling and therefore a large CPE. 
It is a formidable task to figure out in such a situation which G-matrix diagrams should 
be embodied in the Q-box. The two leading terms are the well-known first-order G- matrix 
diagram and the second-order core polarization diagram G^^pih ||TT|]. Hjorth- Jensen et al. |T^ 



have investigated the third-order Q-box diagrams for the sd shell. They concluded that, 
after folding, the net effect on V^ff was a change of about 10 — 15%, as compared with the 
case when only the first- and second- order ones were considered. Higher and higher-order 
core polarization diagrams rapidly become prohibitively more difficult to deal with. Thus 
in practice one can only include some low-order diagrams for the calculation of the Q-box. 
Besides, when there are disagreements between theory and experiment, one is not sure if 
they are due to the NN interaction or to the approximation adopted in solving the many- 
body problems (such as the neglecting of the higher-order core polarization diagrams). The 
environment in the halo nucleons is different and may be more promising, because they are 
located quite far away from the core. A schematic comparison between normal and halo 
nuclei is given in Fig. 1. A relatively weak CPE is expected for halo nucleons, and therefore 
an V^fi predominately governed by the free NN interaction. That is, K// should be in essence 
the bare G-matrix, which presently can be calculated to a high degree of accuracy. Hence the 
halo nuclei, besides having excitingly interesting and exotic properties, may furnish as well 
a much better testing ground for the fundamental NN interactions, than ordinary nuclei. 

Motivated by the above scenario, we present in this letter a microscopic derivation of the 
V^ff for halo nucleons in ^He and ^^Li, starting from the Paris and Bonn NN potentials and 
using a G- matrix folded-diagram approach @,|l3l • The main steps in such a derivation are: 



i) Choice of the model space P. An important criterion for selecting the model space 
P is that its overlap, with the physical states under consideration, should be as large as 
possible. For instance, the ^ife, i.e., the ^He core, should remain essentially as an ordinary 
a-particle, with little perturbation from the distant halo nucleons. For the P space we shall 
use a closed (Osi/2)^ core (a-particle) with the valence (halo) nucleons confined in the Op 



shell. Yet the halo nucleons have a much larger r.m.s. radius than the core, and therefore 
an oscillator constant huj considerable smaller than that given by the empirical formula 
hoj = 4:5 A~^^^ — 25 A~'^^^ MeV (valid for ordinary nuclei). It would not be then feasible to 
reproduce both radii, using shell model wave functions with a common hu. One may get past 
this difficulty by including several major shells in the one-frequency shell model (OFSM) 
calculation. But this would be very tedious. A convenient and physically appealing solution 
to this problem is to employ a two-frequency shell model (TFSM) for the description of halo 
nuclei, as suggested in Ref. p. Within the TFSM one uses oscillator wave functions with 
huin and huJout for the core (inner) and the halo (outer) orbits, respectively. The notations 
bin and bout also will be used from now on, with 6^ = h/muo. In the present work 6j„ is 
fixed at 1.45 fm, while bout is treated as a variation parameter (or generator coordinate). To 
assure the orthonormality, we have actually used 6j„ for all the £ = waves (Osi/2, lsi/2, ■ ■ ■) 
and bout for waves with other I values. 

ii) Evaluation of the model-space G-matrix. For ordinary nuclei, the G-matrix can be 
calculated rather accurately with the method developed in Refs. |T3|,|T3. We extend below 



this method to the halo nucleons in the context of the TFSM. For a general model-space P, 
we define the corresponding Brueckner G-matrix by the integral equation IT^JTH 



G{uj) = V + VQ2 ^——Q2G{uj), 

UJ - Q2TQ2 

where u is an energy variable, Q2 is a two-body Pauli exclusion operator, and T is the two- 
nucleon kinetic energy. Note that our G-matrix has orthogonalized plane-wave intermediate 
states. The exact solution of this G-matrix is G = G^? + AG [[TSlJT^ , where Gp is the free 
G-matrix, and AG is the Pauli correction term 

with e = u — T. The projection operator P2, defined as (1 — Q2), will be discussed later. 
The basic ingredient for calculating the above G-matrix is the matrix elements oi Gp within 
the P2 space. This space contains all the two-particle states that must be excluded from 



the intermediate states in G-matrix calculations. For ordinary nuclei, where the OFSM 
is used, the states excluded by the Pauli operator and those contained within the model 
space have a common length parameter b. For halo nuclei, where we use the TFSM, the 
situation is more complicated, as the wave functions for the excluded states and those 
within the model space have in general different length parameters 6j„ and bout- Hence 
to calculate AG, we need the matrix elements of Gi? in a bm — bout mixed representation. 
This poses a technical difficulty because the transformations, from the cm. coordinates to 
the laboratory coordinates for two-particle states with different oscillator lengths, are not 
as easy to perform as for one common oscillator length. We have adopted an expansion 
procedure to surmount this difficulty. Namely, we expand the oscillator wave functions with 
bin in terms of those with bout, or vice versa. When 6j„ and bout are not too different from 
each other, this procedure is relatively effortless to carry out. Usually a high accuracy can 
be attained by including about 8 terms in the expansion. Still, the calculation of the two- 
frequency G-matrix is significantly more complicated than the ordinary one-frequency one. 
Another difficulty, in deriving the G-matrix for halo nuclei, is the treatment of its Pauli 
exclusion operator. As the halo nucleons are rather far from the core nucleons, the effect of 
Pauli blocking is expected to be small. But, to get a reliable result for a small effect, a very 
accurate procedure has to be employed. We write the projection operator Q2 as 

Q2 = E Qiab)\ab){ab\, 

all ab 

where Q{ab) = 0, if 6 < rii, a < n^, 01 b < n2, a < n2, or b < n^, a < rii, and Q{ab) = 1 
otherwise. The boundary of Q{ab) is specified by the orbital numbers (ni, n2, ^3). We denote 
the shell model orbits by numerals, starting from the bottom of the oscillator well: 1 for 
orbit Osi/2, 2 for 0^3/2, ■ ■ ■ , 7 for O/7/2, and so on. rii and ^2 stand for the highest orbits of 
the closed core (Fermi sea) and of the chosen model space, respectively. For example, we 
consider ^He as a closed core and all 6 orbits in the sp and sd shells are included in the model 
space. Then ni = 1 and n2 = 6. As for the G-matrix intermediate states we consider only 



particle states {i.e., states above the Fermi sea), in principle should be 00 |jT^. Still, in 



practice this is not feasible, and has to be determined by an empirical procedure. Namely, 
we perform calculations with increasing values for until numerical results become stable. 
In Table 1, we display some representative results of our two-frequency G-matrix for the 
{OsOp} model space, with hm = 1.45 fm and bout = 2.0 fm. The only approximation here 
is the finite truncation. A satisfactory convergence is attained for = 21, and this 
value is used here. It is worth noting that, although the halo nucleons are widely separated 
from the closed core, the Pauli correction term AG {= G — Gp) is still quite significant. 

iii) Calculation of the irreducible diagrams for the vertex function Q-boxes. Once derived 
the G-matrix, we can calculate the irreducible vertex function Q-hox. Finally, the model- 
space energy-independent K// is evaluated in terms of the Q-hox folded-diagram series P,p!0 



following closely the procedures of Ref. |11 



Diagonal matrix elements of G, G-^pih and K//, for the states 1(^3/2)^; T = 1, J = 0) and 
\{pi/2)'^',T = 1, J = 0), are shown in Fig. 2 as a function of bout, for both the Paris and 
Bonn-A potentials. As we increase bout, we are augmenting the average distance between 
the halo nucleons and the core and so reducing the coupling between them. For sufficiently 
large bout, the total CPE must be small and it should be sufficiently accurately given by the 
second-order (lowest order) core polarization diagram alone. In fact, as bout increases, the 
core polarization diagrams Gjpih approach rapidly and monotonically to zero and become 
negligibly small at bout — 2.25 fm. In our TFSM approach, we have assumed a fixed ^He 
or ^Li core, always described by 6j„ = 1.45 fm. Therefore the energy denominator for the 
diagram G^pih is fixed by the corresponding core and does not change with bout- This means 
that the suppression of G^pih is entirely due to the weakening of the core-valence particle 
interaction. The behavior of the bare G-matrix and V^ff shown in Fig. 2 are also of interest. 
First, they are quite similar to each other. Second, while for the P3/2 case, they become 
weaker as bout increases, in the pi/2 case they become stronger as bout increases. Third, at 
large bout the results given by the Paris and Bonn A potentials are practically identical. This 
is because their long-range parts do not differ much from each other. 

To assess to which extent the nuclear model formulated above is reliably it is necessary 
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to compare our results with experiments. The valence or pairing interaction energy between 



the halo nucleons in He is obtained from the odd-even mass difference |T6| 



E^^'Pi^He) = -[B^He) + B^He) - 2B{^He)] = -2.77 MeV. 



To calculate this energy, we need to diagonalize the folded-diagram V^fj [0] in a ps/2 and 
Pi/2 model space (T = 1, J = 0). In this way we get E^^ = —2.97 MeV at bout = 2.25 fm 
for the Paris potential. As the CPE is strongly suppressed for such a bout value, this result 
almost entirely comes from the bare G-matrix. The ground-state wave function of ^He is 
almost pure (^3/2)^, with very little admixture. Thus our Ejl^ is also close to the 

diagonal G-matrix element shown in Fig. 2. As we have used a P3/2 —Pi/2 model space, our 
wave function for ^^Li has only one component (neutron orbits closed). Then the diagonal 
(pi/2)^(r = 1, J = 0) matrix element of K//, which is quite close to the unfolded value of 
Fig. 2 at bout = 2.25 fm, is directly comparable to the pairing energy for ^^Li. From the 
masses of ^^Li, and %i [|1^], we get E^^p = —1.14 MeV, while our result is -E*^ = —0.81 
MeV at bout = 2.25 fm for the Paris potential. This discrepancy could be pointing out 
that some physics is still missing in our description of the valence ^^Li neutrons. It is very 
likely that they should not be entirely confined to the p shell, but a larger space, such as 
{0pi/20dls}, is probable needed. 

In passing, we mention that the above bout = 2.25 fm is a reasonable choice. Recall that 
we have fixed b^ = 1.45 fm. With these values of 6j„ and bout, and assuming a pure s'^p"' 
wave function, we get that R^^i^He) = 2.51 fm, in good agreement with the empirical value 
R^'^P^^He) = 2.57 ± 0.1 fm [0; similarly, R^^{^^Li) = 3.03 fm while R'^^'Pi^^Li) = 3.1 ± 0.1 



fm [0- 

We have also calculated the valence interaction energy for ^Li using a similar folded- 
diagram procedure in the P3/2 — P1/2 space. From the empirical masses of ^Li, ^Li, ^He and 
^He |T6[, we obtain E^^p = —6.56 MeV. (This number was incorrectly given as —3.55 MeV 
in Ref. [Q.) Our result is E^^ = —6.64 MeV for the Paris potential, if we use bout = 1-75 
and bin = 1-45 fm. It is of interest to stress that ^Li is not a halo nucleus, according to 



our calculation, in the sense that there is no need to employ a very large bout for its valence 
nucleons. 

In summary, we have derived the effective interaction for the valence nucleons in halo 
nuclei, starting from realistic NN interactions. Our preliminary results are encouraging. We 
have employed a two-frequency shell model approach, to give a good spatial description for 
both the core nucleons and the halo nucleons. While keeping the inner length parameter 
bin fixed, we gauge the spatial extension of the halo nucleons by varying the outer length 
parameter bout- In this way we have explicitly proved that the core polarization effect is 
strongly suppressed at large bout values, as required by the large empirical r.m.s. radii of 
halo nuclei. Ergo, the effective interaction between the halo nucleons is predominantly given 
by the bare G-matrix alone, in accord with our expectations. The Pauli blocking effect on 
the G-matrix has been found to be very important, and it can be calculated quite accurately 
as we have demonstrated. Thus it appears that one can derive the effective interaction for 
halo nuclei much more accurately than for ordinary nuclei. We enthusiastically believe that 
the halo nuclei, which have already greatly enhanced our knowledge about nuclei, may in 
addition provide a more accurate testing ground for the fundamental NN interaction, than 
the ordinary nuclei. 
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TABLES 

TABLE I. Dependence of the two-frequency G-matrix on the choice of ns. Listed are the 

matrix element {{0p^/2)'^;TJ\G{ll!)\{0p3/2)'^iTJ) (in MeV), calculated for the Paris potential and 
three different values of lo (in MeV), with TJ = 01 (upper panel) and with TJ = 10 (lower panel). 
We have used 6,„ = 1.45 and bout = 2.0 fm for the length parameters, and ni = 1 and n2 = 6 for 
the exclusion operator. The first row in each group (F) denotes the free G-matrix. 

n3 Lo = —5 u = —10 Lo = —20 

F -6.896 -4.530 -3.155 

6 -2.218 -2.115 -1.885 

15 -2.217 -2.114 -1.882 

21 -2.217 -2.114 -1.882 



F -4.422 -3.933 -3.480 

6 -2.768 -2.748 -2.701 

15 -2.761 -2.744 -2.698 

21 -2.761 -2.744 -2.698 



Figure Captions 

Fig. 1 Comparison of core polarization in ordinary and halo nuclei. 

Fig. 2 Diagonal matrix elements of G^pih (dotted lines), G (dashed lines) and V^ff (full 
lines) for the states |(p3/2)^;T = 1, J = 0) (upper panel) and |(pi/2)^;T = 1, J = 0) (lower 
panel) as a function of bout] calculations done with Paris and Bonn- A potentials are shown 
by open and solid symbols, respectively. The G-matrix curves are for cu — —5 MeV and 
Pauli exclusion operator with (711,712,^3) = (1,3,21). 
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Halo Nucleus 



